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Abstract 

We consider single particle Schrodinger operators with a gap in the energy spec- 
trum. We construct a countable set of exponentially decaying functions, which form 
a complete, orthonormal basis set for the states below the spectral gap. Each such 
function is localized near a closed surface. Estimates on the exponential decay rate 
and a discussion of the geometry of these surfaces is included. 

1 Introduction 

The fast developing field of nano-science and the need for microscopic understanding of 
the biological processes are just two of the driving forces for the search of computationally 
efficient electronic structure algorithms. One important goal is to develop electronic structure 
algorithms that scale linearly with the number of particles. It seems that such algorithms 
will involve a localized, real-space description of the electronic structure [1] and, of course, 
Wannier functions will play a major role [2, 3]. 

The Wannier functions were originally defined for perfectly periodic crystals [4]. Their 
key property is the exponential localization, which has been proven for ID periodic systems 
[5] and for simple bands in more than 1 dimension [6, 7, 8]. 

A natural and important question is if one can define the equivalent of the Wannier 
functions for systems that are not periodic. Examples of such generalized Wannier functions 
were given by Kohn and Onffroy [9] for crystals with one impurity and Rehr and Kohn [10] 
for crystals with surfaces. In a relatively recent paper [11], Nenciu conjectured that any 
system with a gap in the energy spectrum should posses exponentially localized Wannier 
functions and he proved the conjecture in one dimension. His proof is based on an idea 
introduced by Kivelson [12], who showed that, for ID perfectly periodic systems, the eigen- 
vectors of PqxPq provide a set of Wannier functions (Pq = the band projection operator). 
Moreover, it was later shown by Marzari and Vanderbilt [13] that these Wannier functions 
are maximally localized. Now, if one tries to extend Kivelson's construction to more than one 
dimension, the problem is that, even for perfectly periodic systems, PqxPq, PquPq and PqzPq 
do not commute, and any of these three operators have, in general, continuum spectrum 
and, consequently, no localized eigenvectors. 

In this paper, we propose an alternative way of looking into the problem. Let us first state 
our main result. Consider a V^-bounded potential v over M*^ {d = 1,2, .. . ), with relative 
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bound less than one. We consider only potentials that are bounded from below. Measuring 

the energy from the bottom of the potential, it is equivalent to say that the potentials are 
positive. Other than this, the only assumption on v is that the self-adjoint Hamiltonian, 

: P(-V^) ^ L^(R'^), H^-V^ + v, (1) 

has a gap in the energy spectrum. Let /C be the space of all the states below the gap. Our 
goal is to construct a countable set of orthonormal functions, all decaying exponentially and 
forming a complete basis set for /C. The solution is given in the following theorem. 

Theorem 1. For q > 0, let Wq be the following bounded, self-adjoint operator: 

Wq-.IC^JC, Wq = Poe-'^fPo, (2) 

where p = [1 + r^]^/^. Then the following are true. 

i) The spectrum of Wq, a{Wq), is discrete. 

ii) All eigenvectors decay exponentially (with the same rate) as |r| — > oo. 

Since Wq is self-adjoint and <j{Wq) is discrete, its eigenvectors form a complete and 
orthonormal basis for /C. Thus, simply by solving for the eigenvalues of Wq, one can construct 
a complete (in /C), orthonormal and exponentially decaying set of functions. These functions 
will not have all the localization properties listed in Ref. [11]. Instead we can prove that the 
eigenvector corresponding to the eigenvalue Aj is exponentially localized around the sphere 
of radius 

and centered at the origin. We believe that such an eigenvector is localized around one or 
more points on the sphere and, in general, the localization in the direction normal to the 
sphere is exponential while the localization parallel to the sphere is not. However, these 
functions can be as useful as the usual Wannier functions. For example, the matrix element 
of the Hamiltonian, between a state localized around the sphere of radius Ti-m and a state 
localized around the sphere of radius ri+„, decays exponentially with m and n. 

By simply modifying the function form of p, we can obtain generalized Wannier functions 
that are localized around more general closed surfaces (like ellipsoids, etc.). 



2 Proof of the main statement 

The proof of our main statement involves the analytic continuation technique, first used 
by des Cloizeaux [6, 7] in the context of band calculations. Let Pq denote the projection 
operator onto the states below the gap. Let E± denote the upper respectively the lower 
edge of the gap A and consider the Tinitary transformation e~*^^, with q real. We define the 
following unitarily equivalent representations of the band projector operator: 

Po{q) = e-^^^'Poe'^". (4) 
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We will show later that we can analytically extend this family to complex q. More precisely: 
Lemma 1. Pq (?) extends to an analytic family of bounded operators, 

ll^o(?)|| <oo, (5) 

for any g in a complex strip |Img| < qm = sup ^/ {Aij:, Aij:), for any operator A). 

Moreover, qm is always larger than 



(6) 



where qo = \/E_ + \/E^. 

We start now the proof of Theorem 1. i). We write as 



W,^P,\ e-^'^iz-H)-'^^, (7) 



r 



with r a contour surrounding the lower band. Since e^'^f{p'^ + a)^^ is Hilbcrt-Schmidt [14] 
and (p^ + a){z — H)^^ is bounded, it follows that Wq is Hilbert-Schmidt. As a consequence, 
(T{Wq) is discrete. Moreover, the eigenvalues accumulate at zero and only at zero. 
a) Let be a normahzed eigenvector, Wq-^ — Xt/j (of course, PqV' — '0) ■ Since 



A 



= J e-'i''\iP{f)\'^df (8) 



and p > 1, all eigenvalues are strictly positive and smaller than e^^. 

Let us consider first the case when q is smaller than the qM of Lemma 1. In this situation 
we can prove 



/ 



e2«^|V'(r)|W< oo, (9) 



which clearly shows that \ip{r)\ decays faster than e as |r| — > oo. Indeed, observing that 
eiPWq = Po{iq)Po, we find 



1/2 



= We'lfipW = X-^\\e'ifWqip\\ (10) 
= X-'\\Po(iq)^\\ ^ X-'\\Po{iq)\\, 



and, by definition, qM is the maximum value of q for which ||Po(^?)|| < oo- 
For the case when q ^ qu, one can show that 



J e^''f\iP{f)\^df< oo, (11) 
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for any k < qm- Indeed, 

||e'=''^|| = A"^ ||Po (ik) e-^^-'^)"^!! ^ A"^ ||Po {ik)\\ < oo (12) 

and this completes the proof of Theorem 1. 

Thus, any Wq with q > provide a complete and orthonormal set of exponentially 
decaying functions. The exponential decay rate is larger or equal to the minimum between 
q and qM- 

Wc now discuss their localization. We will limit ourselves to the case q < qM- Let us 
write the eigenvalues of Wg as Aj = e"^''*. Since Aj are positive and accumulate at zero, 
{Pj}^ is a sequence of positive numbers which accumulate at infinity. If ipi is an eigenvector 
corresponding to Aj, then we have from Eq. (8), 

y e''(^^-'')|V'i(r)|'(if = 1 (13) 

and from Eqs. (10), 

j e'«(''-''^)|^i(r)|'(if ^ \\Po{iqW- (14) 

Observing that the right side of the last two equations is independent of the index i, we can 
conclude that ipi is exponentially localized (to annihilate the exponential factors) around the 
sphere defined by p = pi, i.e. the sphere of radius rj given in Eq. (3). More precisely, ipi 
falls off exponentially with a rate q as p — pi ^ oo and with a rate q/2 as p — pj — > — oo. Of 
course, the maximum localization is obtained when q = qM- 

The rest of the paper contains the proof of Lemma 1, which provides a lower bound for 
qM- This lower bound can be easily estimated from the band structure and we found, from 
numerical calculations, that qc is a very good approximation of qM for small gaps. 

3 Proof of Lemma 1 

The proof follows Barbaroux et al [15]. We will try to ehminate the unnecessary constants 
and optimize the technique. 

Proposition. Consider A and -B, two bounded, self-adjoint operators. Assume that A has 
a spectral gap located at 0, 

d± = (iisi((7±(A),0) > 0, (15) 

where (t± denote the parts of the spectrum located above and below zero. Let P± denote 
the spectral projector onto the states corresponding to c±(^) and C > be a bounded 
self-adjoint operator such that s+ = d+ — ||P+CP+|| > 0. Then, for q real, 

\\{A-C + iqB)-'\\^[l-'^^^ 'niax|^,^|, (16) 
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where 



'^jp:bpj\- ^^^^ 



Proof. Let be a norm one but otherwise arbitrary vector, ip± = P±(p and 5+- = Pj^BP^. 
We start from the observation 

Re ((^+ - (A - C + iqB)ip) ^ - C + iqB)^p\\ . (18) 

The left hand side is equal to, 

(<^+, {A - C)99+) - ((^_, {A - C)^_) - 2glm (<^+, B^_) , (19) 

and we have successively: 

\\{A-C + iqB)^\\ (20) 

1 1 2 1 1 2 

= (1 - kl /?) (s+ llv^+f + ||v9_f ) 

+ kl \\B^-\\{0_\WA\-^,y-\\)\ 

> (1 - |g|/g)min{s+,ci_} (21) 

and the affirmation follows. ■ 

We define the following analytic family, 

Hq = e-''iPHe''iP + q^\Vp\^ + g (Vp • p + p- Vp) , (22) 

where we notice that |Vp| ^ 1. Since V p-p + p- Vp is relatively if-bounded, Hq extends to 
an analytic family of type A over the entire complex plane [16]. Then, according to Combes 
and Thomas [17], 

P„(,) = ^(.-i/,)-'^. (23) 

where F surrounds the lower band. One can sec that Po(g) belongs to an analytic family of 
bounded operators as long as V belongs to the resolvent set of Hq. As we will argue later, 
we have to worry only about the point where F crosses the real axis. We consider then an 
energy E in the spectral gap and calculate for which values of q the resolvent {Hq — E)~^ 
remains a bounded operator. It is enough to consider only purely imaginary q. In the above 
Proposition, we take 

A^\H - E\-^/'^{H - E)\H - Ey-^/"^, (24) 
B^\H -E\-^''^{Vp-p + p-Vp)\H -E\-'^l'^ (25) 
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and 



C = q^\H - E\-^'''\V p\''\H - E|-V2. (26) 
Then we have — 1 and s+ = 1 — — E), and we evaluate from, 



15. 



sup 



+ {Vp-p\H- E\-^/^<p+, \H - E\-^l^ip_ 



(27) 



where the supremum goes over all norm one 99+ and ip- in the upper respectively lower band. 
If 



^l:± = {H + af'^\E-H\-^'^^^, 
for some positive a, then we find 



< \\p{H + ay 



-1/21 



+ 



+1 



•\J E^ — E ^ E — E_ j 



(28) 



(29) 



Since the potential is positive, ||p(-ff + a) ^''^H < 1 for all a > 0. We take the limit a = 
and use the spectral theorem to evaluate HV'itll and finally find 



\B^-\\ < 



^ + y/E: 



^{E+-E) {E-E_y 
From Eq. (16), we can conclude 

\H - E\^'^ {Hiq - Ey^ \H - E\^^^ 



where 



s+(i - \q\/q)'' 

y/{E-E_){E+-E-q^) 
JEZ^JEl 



(30) 



(31) 



(32) 



The optimal position of the energy E in the spectral gap for q to become maximum is 
E^ = (£;_ + £;+ - g2)/2 which leads to g = (A - (f')l1qQ. Then from Eq. (31), 



and the right side is finite as long as 



q<qA x/l + A/go'-l . 



(33) 



(34) 



As one can easily see, adding an imaginary part to Ef. will not affect the real part in 
Eq. (18) and the immediate consequence of this is that {Hig — Ec + ie)~^ is bounded as 
long as {Hiq — Ec)~^ is bounded. Consequently, for |lmg| < qc, we can always find a curve F 
surrounding the lower band and lying in the resolvent set of HiqM 
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